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ABSTRACT 



Doppler observations from Keck Observatory have revealed a triple planet 
system orbiting the nearby mid- type K dwarf, HIP 57274. The inner planet, 
HIP 57274b, is a super-Earth with Msini = 11.6 Me(0.036 Mj^p), an orbital 
period of 8.135 ± 0.004 d, and slightly eccentric orbit e = 0.19=b0.1. We calculate 
a transit probability of 6.5% for the inner planet. The second planet has Msini 
= 0.4 Mjup with an orbital period of 32.0 ±0.02 d in a nearly circular orbit, and 
e = 0.05 ±0.03. The third planet has Msini = 0.53 Mjup with an orbital period 
of 432 ±8 d (1.18 years) and an eccentricity e = 0.23 ± 0.03. This discovery 
adds to the number of super Earth mass planets with Msini < 12M® that 
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have been detected with Doppler surveys. We find that 56 ±18% super-Earths 
are members of multi-planet systems. This is certainly a lower limit because 
of observational detectability limits, yet significantly higher than the fraction of 
Jupiter mass exoplanets, 20 ± 8%, that are members of Doppler-detected, multi- 
planet systems. 

Subject headings: planetary systems - stars: individual (HIP 57274) 



Introduction 



Low mass K and M dwarf stars are important targets for exoplanet surveys because of 
their proximity and prevalence in the Galaxy. Differences in the number and type of exo- 
planets orbiting these stars relative to more massive stars reflect conditions in the protoplan- 
etary disk that are important for planet formation theory. Microlensing surveys suggest that 



both ice and gas giant planets are common at separations beyond the ice line ( IGould et al. 



20101 ). However, the fraction of gas giant planets det ected inside the ice line by Dopple r 
surveys is relatively low for late K and early M dwarfs (lEndl et al.ll2003l : iButler et al.ll2006l ). 
Fischer et al.l (120111 ) find that the population of giant plan ets has a precipitous d ecline for 
stars redward of B — V= 1.1, a spectral type of about K5V. ICumming et al.l (120081 ) estimate 
that relative to FGK stars, M dwarfs are far less l ikely to harbor gas giant planets with 
periods shorter than 5 years. Ijohnson et al.l (2010) find that 3.4lQ'g% of low mass stars 



(M < O.6M0) have planets with Msini > O.SMjup and semi-major axes less than 2.5 AU. 

The remarkable discovery of more than 1200 planet candidates by the Kepler Mission 
(IBorucki et al.ll201ll ) provides statistics for smaller planets, and suggests that the reduced 
planet occurrence with later spectral type only applies to gas giant planets. After cor- 
recting for th e poore r detectability of transits around higher mass stars with larger radii, 
Howard et al.l (120111 ) find that 20 - 30% of low mass stars have planet candidates with 
Neptune-like radii between 2-4 while the fraction of more easily detected Jupiter-radii 
planets hovers at a few percent. Howard et al. also see evidenc e of a rising occurrence 
of sm all-radius planets among cooler, less massive stars. Further, ISchlaufman &: Laughlin 
(I2OIII ) find that while the planet-metallicity correlation among sun-like stars is strongest for 
those hosting large-radius planets, the planet-metallicity correlation among low-mass stars 
is significant even among hosts of small-radius planets. Doppler surveys of nearby, low-mass 
stars provide a means of testing whether these correlations hold among stars in the solar 
neighborhood, and if so, inform the target lists of future planet search efforts. 



Late type stars are especially appealing targets for rocky planet searches in Doppler 
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surveys because the lower stellar mass results in a larger reflex velocity for a given mass 
planet. Fu rthermore, chromospheri c activity in low rn ass stars has less impact on the radial 
velocities (llsaacson &: Fischerll2010l : iLovis et al.ll201l[ ). The ubiquity of low mass stars, cou- 
pled with more easily detected Doppler signals and lower stellar jitter all make K and early 
M dwarfs desirable targets in the search for rocky planets. However, a caveat has emerged: 



the in ner planetary architectures of low mass stars may be more complex. iLatham et al. 



( I2OIII ) analyzed multi-planet systems detected in transit with the Kepler Mission and find 
that solar type and hotter stars are more common hosts of single transiting planets, while 
multi-planet systems are more often detected around cooler stars. Among the 170 multi- 
planet systems detected by the Kepler Mission, 78% contain planets no larger than Neptune ; 
close-in gas giant planets are far less common in multi-planet systems (jhatham et al. 2011 ). 
This has profound implications for Doppler surveys: the challenge of detecting the small 
velocity amplitudes of Neptune-like planets will be compounded by the need to deconvolve 
multiple signals of similar amplitude. For both of these reasons, a larger number of observa- 
tions over a longer interval of time are required to resolve the components of these planetary 
systems. 



To better understand the 



requency and ar chitectures of planetary systems around low 



mass stars, we began "M2K" (lApps et al.ll2010l ). a Doppler survey of M a r id K dwarf stars 



draw n from the SUPERBLINK proper motion survey Jlepine fc Sharalboosi iLepine fc Gaidod 
201lh . Here, we report the detection of a triple planet system orbiting HIP 57274 comprised 
of a super Earth-mass planet and two planets that are likely gas giants. 



Observations and Data 



Doppler observations are carried out with the Keck 10-m telescope and the HIRES 
spectrograph (IVogt et al.l 119941 ). An iodine cell is used to provide the wavelength solution 
and sampling of th e line spread function to model the Doppler shift in the stellar spectra 
( iButler et al.lll996[ ). The B5 decker on HIRES provides a spectral resolution of about 55,000 
and an exposure meter terminates the observations when a target signal-to-noise ratio of 
about 200 is achieved. Most of the M2K stars are fainter than = 9, requiring exposure 
times of 10 - 15 minutes. We have acquired three or more Doppler measurements for more 
than 170 stars, with formal measurement uncertainties of about 1.5 m s~^. 



-4- 



2.1. HIP57274 



HIP 57274 (GJ 439) has an apparent magnitude of V 



(ESA 


1997; 


van Leeuwen 


2007) 



This yields a distance of 25.9 pc and absolute visual magnitude of My = 6.89. We carried out 
spectr al synthesis modeliiig of t h e iodine-free template sp ectrum using Spectroscopy Made 
Easy (jValenti fc Piskunovl Il996l : IValenti fc Fischer to determine stellar parameters. 



Following the method described in IValenti et al. 



(120091 ) ■ the initial param eters derived with 



SME were used as input for interpolation of the Yonsei-Yale (Y2) isochones (jPemarque. et al. 



20041 ). which returns a new value for log^r. We then ran an iterative loop, fixing logs' to the 
isochrone value and running a new SME model fit. The other (free) stellar parameters change 
in response to the fixed surface gravity, so subsequent isochrone interpolations produce a 
slightly different value for logg. We continue the iteration until the output logs' form the 
isochrones does not change by more than 0.001 dex from the previous iteration. This provided 
the following spectroscopic parameters: T^s = 4640 ±100K, v sinz = 0.5 ±0.5 km s~^, [Fe/H] 
= +0.09 ±0.05, logs = 4.71 ±0.1, and Y2 isochrone models for a stellar mass of 0.73 ±0.05 
Mq, radius of 0.68 ±0.03 Rq and an age of 7.87 ±5 Gyr. The spectral classification is listed 
as K8V in the Hipparcos catalog, although the B — V color, spectroscopi c temperature and 
derived mass are more consistent with a slightly earlier spectral type. I Gray et al.l (120031 ) 
list a spectral classification of K4V for this star. A medium-resolution spectrum, obtained 
by one of us (SL) with the Mark III spectrograph at the MDM 1.3-m telescope, returns a 
spectral type of K5V, and we adopt this spectral classification for this paper. The stellar 
parameters for HIP 57274 are summarized in Table 1. 



2.2. Chromospheric Activity and Velocity Jitter 



Isaacson fc Fisched (12010! ) determined the chromospheric activity of 2630 stars observed 
at Keck by measuring the emission in the cores of the Ca II H & K lines relative to adjacent 
continuum regions. These 5'hk values were cal ibrated to the long-standing S'hk values from 
the Mt. Wilson program (jPuncan et al.lll99ll ). Using their technique, we measure a mean 
•S'hk = 0.38 for HIP 57274. The individual measurements of S'hk are listed in the last column 
of Table 2, along with the radial velocity measurements. 

Cooler stars typically have larger S'hk values than solar type stars because of weaker 
continuum in the near-UV. Therefore, S ' hk va lues should not be directly compared for stars 
of different spectral types. iNoyes et al.l (119841 ) correct for the photospheric contributions to 
produce a normalized activity metric, logi?^j^. Chromospher ic activity is tied t o dynamo- 
driven flows and decreases as the star ages and spins down. iNoyes et al.l (119841 ) made use 
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of this activity-rotation correlation and calibrated log i?^K to rotational periods for stars in 
open clusters of different ages. Using their relation, we derive log = -4.89, indicating low 
activity for HIP 57274 and Prot ~ 45 d, consistent with the SME-estimat e d age of ~ 8 Gyr. 
We caution that both log ii^K Prot were calibrated by iNoyes et al.l ( 119841 ) using stars 
with 0.4 < B — V < 1.0 and rotational periods shorter than 30 days; HIP 57274 falls outside 
both of these properly calibrated ranges and therefore our derived logi?^^: rotational 
period should be considered to be less certain estimates. 

Because the activity calibration for stars redward of B — V = 1.0 is an extrapolation, 
Isaacson Sz Fischer J2010h established a differential activity measurement, AShk, and eval- 
uated the impact of chromospheric activity on radial velocities in four separate ranges of 
B — V. Following their method, we plot the 6'hk index for the 170 stars observed on the 
M2K program with B — V color from 0.8 to 1.6 (Figure 1). The dashed line in this plot is 
taken from taken from llsaacson fc Fischen (120101 ) and indicates the baseline Suk values for 
low activity stars. AShk is the difference between this baseline activity level and the mean 
S'hk for a given star. Active stars float high above the baseline values while chromospherically 
quiet stars are closer to the dashed line in Figure 1. 

In Figure 2, we plot velocity rms as a function of excess chromospheric activity, AShk 
for all of the observed M2K stars. We fit a linear function to the lower twentieth percentile 
velocity scatter and interpret this (red solid line overplotted in Figure 2) as the quadrature 
sum of internal errors and jitter (where jitter is both instrumental and astrophysical) . The 
rms scatter for inactive stars with AShk ~ 0.0 is 2.38 m s~^, and given the typical internal 
errors of 1.5 m s~^, this implies a minimum jitter of 1.45 m s~^. We measure AShk = 0.03 
for HIP 57274, suggesting a low stellar jitter of ~ 1.5 m s~^. 



2.3. Doppler Observations and Keplerian Model 



Figure 3 shows the time series data, overplotted with our Keplerian mode 



erian mode 
algorithm ( 


for a triple 


Wrieht & Howard 


(Giguere et al. 


2011 


) to 



model the data. The best fit Keplerian model for three planets includes a trend of -0.026 
m per day. Parame ter uncertai nties were calculated with a bootstrap Monte Carlo anal- 



ysis (IMarcy et al.ll2005l ) in KFME (IGiguere et al.ll201ll ). which iteratively fits the data with 



a best-fit model, then adds the scrambled residuals back to the theoretical velocities before 
refitting. 



The planet wiht the shortest period completes one orbit in 8.135 ± 0.004 days and 
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induces a velocity amplitude of 4.64 ±0.46 m s~^. With a stellar mass of 0.73 M©, we derive 
a planet mass Msinz = 11.6 M^and semi-major axis of 0.07 AU. The orbital eccentricity 
is 0.187 ± 0.10 and the argument of periastron passage u ~ 82°. Because the velocity 
amplitude is small compared to the uncertainties and stellar jitter, the eccentricity for this 
planet is poorly constrained, however the signal itself is unambiguous. Figure 4 shows the 
periodogram of the residual velocities of HIP 57274 after removing the other two planets 
and linear trend described below. We carried out a Monte Carlo test to determine the false 
alarm probability (FAP) of the periodogram power. In this test, 10,000 trials were carried 
out where the best fit (triple Keplerian and trend) model was subtracted and the residual 
velocities were scrambled before being added back to the theoretical velocities and refitting. 
In these 10,000 trials, a peak of the same height was never found in the residuals to the 
fit of the two more massive planets, yielding a FAP less than 10~^. Figure 5 shows the 
phase-folded velocities of HIP 57274b overplotted with the theoretical Keplerian model after 
removing the signals from the two additional planets and the linear trend. 

We calcul ated the prospective time of transit, transit duration and transit probability 



using KFME fiGiguere et al.ll201lf ). The transit ephemeris is 2455801.779 ±0.27 HJD, or 
2011/08/28 06:41:40.7 UT. The next transit observable from Mauna Kea occurs at 3 AM 
HST on 13 January 2012, except that the uncertainty in the transit time is more than 6 
hours. The duration of the prospective transit would be 3.08 ±0.356 hours and the geometric 
probability that this planet will transit is 6.5%. 

The middle planet in this system has an orbital period of 32.0 ± 0.02 days. The best 
fit model for this planet has a nearly circular orbit, with eccentricity of 0.05 ± 0.02. The 
velocity amplitude is 32.4 ± 0.6 m s~^, and implies a planet with Msinz = 130 or 0.41 
Mjup and orbital radius of 0.18 AU. The phase-folded data and model for HIP 57274c is 
shown in Figure 6 after subtracting theoretical velocities for the linear trend and the inner 
and outer planets. The prospective ephemeris time is 2454793.035 ±0.176 HJD, although 
the longer orbital period for the middle planet results in a lower 2.7% transit probability. 

The third planet, HIP 57274d, has an orbital period of 432 ±8 days, a velocity amplitude 
of 18.2 ±0.8 m s~^, and orbital eccentricity of 0.27 ±0.03. The planet mass is Msinz =0.53 
Mjup, and the semi-major axis of the orbit is a familiar 1.01 AU. Figure 7 shows the time 
series velocities and Keplerian model for this planet after removing the inner two planets 
and the linear trend. The velocity rms to the triple-planet fit is 3.15 m s~^. However, if we 
adopt the predicted jitter of 1.5 m s~\ we find that = 2.7, indicating that the model does 
not fully describe our observations. 



To check for additional planets in the system, we subtracted the theoretical velocities for 
the linear trend and three Keplerian models. Figure 8 shows a periodogram of the residual 
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velocities, with significant peaks at 1.019 d and 52.996 d that could be ahases of each other, 
due to the diurnal cadence: 1.0 + 1./52.996 = 1.019 d. This predicts the presence of a second 
peak at 1.0 - 1./52.996 - 0.981, and we see a second peak in the periodogram at 0.981 d. 
Interestingly, the ~ 53 day peak has been increasing in strength. Unlike the periodograms 
for the three planets that we modeled (HIP 57472b, c, and d), the 53-d peak is flanked by two 
additional peaks. We suspect that this signal may be caused by spots rotating on the surface 
of the star which reduce flux on the approaching blue-shifted edge of the star and then on the 
receding red-shifted edge of the star. Physically, this would produce a line profile asymmetry 
that could be spuriously modeled as a Doppler shift. We checked to see if the residual 
velocities to the full triple planet fit were correlated with the activity measurements but 
found only an insignificant trend Figure 9 (dashed line). We tried detrending the velocities 
with the linear fit shown in Figure 9, however this only slightly reduced the periodogram 
power in the residual velocities. If we blindly fit this periodic signal with a Keplerian model, 
we derive a period of 53.2 days with an amplitude of 2.6 m s~^, and the residuals drop to 2.63 
m s~^with a xl = 1-28. Interpreting this signal as an additional noise source from coherent 
spots, we can add the 2.63 m s~^ signal in quadrature with the expected jitter of 1.5 m s~^ 
to obtain a revised jitter estimate of 2.8 m s^^. This changes xl to 1.06 for the model of a 
triple planet system plus a linear trend. If the 53 d signal originates from star spots, then 
it should have a detectable photometric signal. We have started a photometric campaign to 
check for spot modulation and to search for the transit signal of the inner 8.135 d planet. In 
addition to the photometric observations, we will continue to obtain Doppler measurements 
to better understand the origin of the 53 day signal. 

We have started a photometric campaign to check for spot modulation and to search for 
the transit signal of the short-period planet. If the 53 d signal originates from star spots, it 
should have a detectable photometric signal. In addition to the photometry, we will continue 
to obtain Doppler measurements to better understand the origin of the 53 day signal. 



2.4. Dynamical stability 

To assess the stability of the HIP 57274 system, we ran an ensemble of dynamical 
sir nulations of the triple planet system with the hybrid symplectic integrator code Mercury 
6 ( IChamberd Il999l ) . Orbital parameters for each body were calculated from the values in 
Table 1 assuming Gaussian-distributed errors (with a truncation at 0.0 for eccentricity). We 
ran 20 simulations with sini = 1 (minimum mass case) and 15 simulations with sini = 0.3 
(95% of random orientations will have a sin i greater than this value) and a mass of 0.73 Mq 
for the central star. For radius and thus collision probabilities, we assumed mean density 
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of 6.0, 1.0, and 1.0 g cm~^ for the b, c, and d planets, which would be typical for the 
densities of rocky and gas giant planets. The time step was set to 0.365 d (4.5% of the 
orbital period of HIP 57274b) and each simulation was run for 10 Myr. In none of the 35 
independent simulations did collisions or ejections occur. We conclude that the system is 
stable regardless of its inclination. 



3. Summary and Discussion 

Here we present a triple planet system orbiting the late K dwarf star, HIP 57274. The 
inner planet orbits in 8.135 days and has a mass, Msinz = 11.6 M^. The orbit is slightly 
eccentric with periastron directed toward our line of site. We calculate a transit probability of 
6.5% with a putative ephemeris times, Tc = 2455801. 776±0. 338 and a duration of 3.08±0.35 
hours. The second planet orbits in 32 days and has Msini = 130.0 M®. The third planet 
has an orbital period of 432 days with a semi-major axis of 1.01 AU and Msini = 167.4 
Mq). The nominal habitable zone of this K star, corresponding to 0.95-1.3 AU around the 
Sun (Kasting et al. 1993) and adjusting for the lower luminosity, lies between 0.41-0.57 AU, 
or orbital periods of 110-180d. No significant periodic signals lie within this range (Figure 
8). 

With the addition of HIP 57274b, there are now 25 planets wit h Msini < 12.0 M^ 



listed in the Exoplanet Orbit Database, or EOD (jWright et al.l 120111 ). It is probable that 
these 25 planets are super-Earths or Neptunes, rather than gas giant planets. Importantly, 
only 8 of the 25 (currently) appear to reside in single planet systems. The remaining 17 
low mass planets are constituents of 10 multi-planet systems. Counting systems instead of 
planets and applying Poisson error bars (i.e., the percentage of single or multiple planet 
systems divided by the square root of the number of planets), we find that 44 ±16% of 
these low mass planets are single (at the current level of Doppler detectability), while 55.6 
±17.6% are in multi-planet systems. To restate: more than half of the Doppler detected 
super-Earths are detected as members of multi-planet systems. 

To compare the architectures of planetary systems containing super-Earths with those 
containing gas giant planets, we extracted all 36 planets from Doppler surveys with M sin i 
between 1.0 - 1.5 Mjup orbiting main sequence stars in the EOD. In this Jovian-mass subsam- 
ple, there were 30 single planet systems (86 ±16%) and 5 multi-planet systems (14 ±6%). 
This result is not particularly sensitive to the arbitrary range of exoplanet mass: in a subsam- 
ple of 63 planets with Msinz from 1.5 - 2.5 Mjup, 73 ±11% systems were single. Another 
sample cut of 41 planets with Msini from 3.0 - 6.0 Mjup yielded 83 ±14% single planet 
systems. Taking an average of these three subsamples, roughly 80 ±10% of the Doppler- 



- 9 - 



detected Jupiters are currently in single planet systems, and about 20 ±8% of Jupit e rs are 
in multi-planet systems. This can be compared with the estimate of IWright et al.l (120091 ) 
who find that at least 28% of planets are in multiple systems. Since they count planets of 
all masses detected before 2009, that result is not inconsistent with our estimate. Relative 
to super-Earths , Doppler surveys detect significantly fewer multi-planet systems with gas 
giant planets. 

In Figure 10, we consider the sibling planets in these multi-planet architectures and 
compare the super Earth and Jupiter subsamples. A total of 26 planets accompany the 10 
super-Earths in known multi-planet systems; these sibling planets also tend to be system- 
atically low in mass - the median Msini = 23.5 M^. In contrast, the Jovian subsample 
includes six companions, spanning a range from 0.58 Mj^p to about 4 Mj^p with a median 
Msini of 2.6 Mj^p. 

Is the dramatic difference in the architecture of low mass planets a bias in Doppler 
detection efficiency or the result of nature (eg., conditions in the protoplanetary disk or 
evolutionary processes)? In Figure 11, we compare the stellar hosts for the super Earth and 
Jovian subsamples. The histogram of stellar mass for the hosts of super-Earths is offset to 
lower mass, with an average of 0.7 Mq, while the host stars of the Jovian sample have a 
mean stellar mass of 1.06 Mq. The dependence of reflex velocity on stellar mass implies 
that the Doppler signals for a star with a mass of 0.7 Mq would be amplifled by about 30% 
relative to a star of 1.06 M^ Therefore, the paucity of Jupiter like planets around lower 



mass stars (IFischer et al.l 1201 ll ) cannot be a selection effect 



However, assessing the scarcity of super-Earths around the more massive host stars of 
Jupiter mass planets is complicated by observational detectability issues. The mean velocity 
amplitude of the super Earth sample would drop from 4.2 m s~^ for the mean host star 
mass of 0.7 Mq to 3.2 m s~^ around solar mass stars. At the same time, as the stellar mass 
i ncreases from 0.7 Mq to 1.06 Mq, the min imum stellar jitter increases by a factor of two 
(llsaacson fc Fischerll2010l : iLovis et al.ll201ll ) or more for chromospherically active stars. As 
a result, if stars with close-in Jovian mass planets also host a system of super-Earths, the 
Doppler signal would be roughly a la detection. Furthermore, the observational biases that 
influence second-planet detection are complex. In some cases the presence of one planet can 
complicate the de tection of additional, lower amplitude planets (for instance if the planets 
are in resonance (lAnglada-Escude et al.ll2010l ) if one planet has poorly constrained orbital 
parameters, or if the observational cadence causes aliased signals near the orbital frequency 
of the additional planet). On the other hand, the presence of a planet can also cause a star 
to receive additional observations in preparation for publication, making the detection of a 
low-amplitude planets more likely. 
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? find that although the migration of giant planes does not completely impede terrestrial 
planet growth, the final accretion phase of terrestrial planets is affected by gravitational 
perturbations from gas giant planets. Although the Doppler detections may only weakly 
constrain the presence of low mass plane t siblings to Jupiter mass planets, the Kepler data 
provide additional support for this case. iLatham et al.l ( 1201 ll ) proposed that the transiting 
Jovian planets detected by Kepler would have migrated into their current locations and likely 
destabilize the orbits of smaller planets. Likewise, nearly all of the Jovian planets detected 
by Doppler surveys are migrated planets that may have driven out close-in Neptunes. 
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Fig. 1. — The emission in the core of the Ca II H & K hnes is a signature of chromospheric 
activity that is parameterized as a S'hk value for each star and plotted ab ove as a function of 
B-V. The dashed line represents the baseline for low activity stars from llsaacson fc Fischer 
(I2OIOI ). Inactive stars fall near the dashed line with AShk ~ 0.0 and active stars have a 
large AS'hk- A filled diamond is used to show the Suk measurement for HIP 57274. 
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Fig. 2. — The rms velocity scatter is plotted as a function of activity. No trends in the velocity 
data were removed. The solid red line is a linear fit to the bottom twentieth percentile velocity 
scatter. Stars with more than 3a jitter (the dashed line) are indicated with green dots and 
represent prospective planet candidates. HIP 57274 is represented by the filled diamond. 
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Fig. 3. — The time series data for HIP 57274 are shown with a Keplerian model for three 
planetary signals plus a linear trend. The model fit has an rms in the residual velocities of 
3.15 m s~^ and xt= 1-06, with an assumed jitter of 2.8 m s~^ added in quadrature to the 
internal errors. 
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Fig. 4. — After subtracting the linear trend and the best fit theoretical velocities for the two 
outer planets, the periodogram of the residuals shows signifiant power at 8.135 days. The 
signal of this inner planet, which we designate as HIP 57274b, has a FAP < 0.0001. The 
two peaks just below and above a period of 1.0 day are aliases of the 8.14 d peak 
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Fig. 5. — The phase-folded radial velocities for HIP 57274b are shown after removing Keple- 
rian signals from the outer two planets and subtracting a linear trend. The Keplerian model 
is plotted with a solid line and has an orbital period of 8.135 ± 0.005 d, orbital eccentricity 
e = 0.19 ± 0.1 and Msmi= 11.6 Mr.. 
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Orbital Phase 

Fig. 6. — The phase-folded radial velocities for HIP 57274c are shown with theoretical 
Keplerian velocities for the inner and outer planets and the linear trend subtracted. The 
Keplerian model has a period of 32.0 ± 0.05 d, eccentricity of 0.05 ± 0.03 and Msinz= 130 
Me or 0.41 Mj,p. 



-20- 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Orbital Phase 



Fig. 7. — The phase-folded data for HIP 57274d is shown with the two inner planets and a 
linear trend removed. The best fit Keplerian model has a period of 431.7 ± 8.5 d, Msinz= 
0.53 Mjup and an eccentricity of 0.27±0.05. 
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Fig. 8. — After subtracting the linear trend and the best fit theoretical velocities for three 
planets, the periodogram of the residuals shows significant power near 1.0193 and 53 days. 
These peaks are likely aliases of each other from the diurnal cadence, and the peak near 53 
days could be associated with spots on the surface of this cool star. 
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Fig. 9. — The residual velocities to the triple planet model are plotted as a function of .Shk 
activity measurements and fit with a first order polynomial and do not show a significant 
trend. 
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Fig. 10. — Excluding the "primary" planets in our two subsamples, we show the distribution 
of sibling planets in the multi-planet systems with super-Earths (red, horizontal hashed 
histogram) and jovian planets (black, diagonal hashed histogram). 
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Fig. 11. — Stellar hosts of detected super Earth planets (red, horizontal histogram) are 
systematically lower in mass than the stellar hosts of Jupiter mass planets (black diagonal 
lines) . 



-25 - 



Table 1. Stellar Parameters: HIP 57274 



Parameter 



V 


8.96 


B-V 


1.11 


Spec Type 


K5V 


Distance [pc] 


25.92 


My 


6.89 


Bol Corr 


-0.50 


L* [Lq] 


0.19 (0.01) 


M, [Mo] 


0.73 (0.05) 


R* [Rq] 


0.68 (0.03) 


Age [Gyr] 


7.87 (5) 


Teff[K] 


4640 (100) 




4.71 (0.1) 


[Fe/H] 


+0.09 (0.05) 


V sin z [km s"''"] 


0.5 (0.5) 


'S'hk 


0.39 (0.02) 




-4.89 


Prot 


45d 


Radial Velocity [km s~^] 


-4.7 
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Table 2. Radial Velocities for HIP 57274 





RV 




Shk 


JD-2440000 


(m s"-*-) 


(m s"^) 




14806.13297 


0.65 


1.34 


0.469 


14807.14718 


-2.97 


1.88 


0.469 


14809.15850 


-1.80 


1.26 


0.468 


15190.15800 


-14.89 


1.18 


0.409 


15198.16272 


-7.39 


1.29 


0.403 


15232.04776 


11.50 


1.26 


0.411 


15255.86086 


-7.77 


1.12 


0.417 


15255.96257 


-7.31 


1.28 


0.415 
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Table 3. Orbital Parameters for HIP 57274 



Parameter 


D 


c 


a 


P(d) 


8.1352 (0.004) 


32.03 (0.02) 


431.7 (8.5) 


Tp - 2.44x10^ (JD) 


14801.015 (1.3) 


15785.208 (9.5) 


15108.116 (14) 


ecc 


0.187 (0.10) 


0.05 (0.02) 


0.27 (0.05) 


u (deg) 


81 (59) 


356.2 (120.0) 


187.2 (5) 


Ki (m s-^) 


4.64 (0.47) 


32.4 (0.6) 


18.2 (0.5) 


Msini (Me) 


11.6 (1.3) 


130 (3) 


167.4 (8) 


ttrel (AU) 


0.07 


0.178 


1.01 


T, (HJD - 2.44e6) 


15801.779 (0.271) 


15793.035 (0.176) 




Tduration (hr) 


3.08 (0.35) 


5.88 (0.1) 




tprob 


6.5% 


2.7% 




Trend (m s'^-^) 


-0.026 (0.002) 






Avg SNR 


225 






RMS (m s-i) 


3.15 






Nobs 


99 








1.06 






Jitter [m s~^] 


2.8 







